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AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 



Length 
Time.. 
Force- 
Power. 
Speed- 



Symbol 



P 
V 



Metric 



Unit 



meter. _ _ 

second! 

weight of 1 kilogram . 



horsepower, (metric) . 
("kilometers per hour_ 
\meters per second - _ 



Abbrevia- 
tion 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour) — 
weight of 1 pound 

horsepower , — 

miles per hour. _ _ 
feet per second — 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or hr.) 
lb. 



hp.:, 

m.p.h. 
f.p.s. 



2. GENERAL SYMBOLS 



W, 
9, 

m, 
h 



Weight==mp 

Standard acceleration of gravity =9. 80665 
m/s 2 or 32.1740 ft./sec. 2 
W 

Mass=; — 
g 

Moment of inertia =mP. (Indicate _ axis of 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



v, Kinematic viscosity ^ ■.'«'■ 

p, Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m _ *-s 2 

15° C. and 760 mm; or 0.002378 lb.-ft."* sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 

0.07651 lb./cu. ft. 



at 



or 



3. AERODYNAMIC SYMBOLS 



s, 

h, 

c, 
p 

s' 
v, 

z> 

L, 
D, 

Do, 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 

Dynamic pressure =^pV 2 

■ L 

Lift, absolute coefficient G L —-^ 



Drag, absolute coefficient G D — ~g 
Profile drag, absolute coefficient ^d 0 =^ 
Induced drag, absolute coefficient @Di—~g 
Parasite drag, absolute coefficient &b P —jj^i 
Cross-wind force, absolute coefficient C c —- 
Resultant force 



Q, 

VI 



o 

qS 



OCf, 
da, 

y> 



Angle of setting of wings (relative to thrust 
line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000 ; or for a model 
of 10 cm chord, 40 m.p.s., the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge, to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
Flight-path angle 
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SUMMARY 

The results of previous reports dealing with airfoil sec- 
tion characteristics and span load distribution data are 
coordinated into a method for determining the air forces 
and their distribution on airplane wings. Formulas are 
given from which the resultant force distribution may be 
combined to find the wing aerodynamic center and pitching 
moment. The force distribution may also be resolved to 
determine the distribution of chord and beam components. 
The forces are resolved in such a manner that it is unneces- 
sary to take the induced drag into account. 

An illustration of the method is given for a monoplane 
and a biplane for the conditions of steady flight and a 
sharp-edge gust. The force determination is completed by 
outlining a procedure for finding the distribution of load 
along the chord of airfoil sections. 

INTRODUCTION 

This report originated in a request of the Bureau of 
Air Commerce, Department of Commerce, for a coor- 
dinated system of applying airfoil section data to the 
determination of wing forces and their distribution. 

The system presented herein yields, within the limi- 
tations of our present knowledge of aerodynamics, a 
general solution of the resultant wing forces and mo- 
ments and their distribution. For the sake of complete- 
ness and facility in use, the report contains a table of the 
important section parameters for many commonly used 
sections and all other necessary data required to solve 
the most practical design problems coming within the 
scope of the system. 

Although the usefulness of the system extends into 
several phases of airplane design, its application to 
structural design is illustrated by following through a 
wing loading condition corresponding to that specified 
in reference 1. 

Two basic principles underlie the system employed. 
First, a force coefficient is treated as the independent 
variable, thus eliminating, as far as possible, the angle 
of attack; and second, the forces are derived throughout 
in terms of certain basic parameters of the airfoil sec- 
tion, which are tabulated for each airfoil section. The 
method followed then builds up the forces progressively 
from simple combinations of certain basic forces and 
simple formulas involving the basic airfoil section pa- 



rameters. As the forces are thus built up, they are re- 
solved into any convenient components. This method 
also has another important advantage in that the in- 
duced drag, which is really only a component of the local 
lift at each section, may be entirely eliminated from the 
analysis. 

In some problems it is desirable to know the location 
of the aerodynamic center of the wing and the pitching- 
moment coefficient about this center in order to con- 
struct the balance diagram of the complete airplane. 
Methods are therefore given for determining these two 
properties. For problems in which the aerodynamic 
center and the pitching moment are not required, a direct 
solution of the forces and force distribution can be 
made. 

BASIC CONSIDERATIONS 

The forces on a wing may be considered to be func- 
tions of the characteristics of the airfoil sections and of 
the spanwise distribution of lift. At a given section 
lift coefficient, the resultant air force and moment on 
the section are, according to wing theory; assumed to be 
independent of all geometric properties of the wing 
except the section shape; moreover, the forces and 
moments acting on any individual section may be 
considered to be independent of adjacent sections or of 
other characteristics of the wing, except as they affect 
the lift distribution and thus the local lift coefficient 
at that section. 

The problem is thus divided into two parts: First, 
the determination of the spanwise lift distribution; 
and, second, the determination of the corresponding 
forces and moments at each section and the summation 
of these quantities to obtain the corresponding forces 
and moments for the entire wing. The spanwise lift 
distribution is obtained in terms of values of the local 
section lift coefficient Ci 0 for a number of sections dis- 
tributed along the span. The subscript zero is used to 
distinguish this section lift coefficient, perpendicular to 
the local relative wind at the section, from the lift 
coefficient c t perpendicular to the relative wind at a 
great distance from the wing. The lower-case letters 
used for these coefficients have been chosen to dis- 
tinguish the lift coefficient for a section (c^dL/qcdy) 
from the usual lift coefficient for the wing, C L . 
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In order to permit easy reference, the symbols used I 
in the text, the figures, and the tables are grouped in 
appendix C. 

For many purposes, it is convenient to express the 
air forces in terms of components along two axes fixed 
with respect to the airplane rather than as the usual 
components, lift and drag. This resolution is con- 
veniently accomplished from the c tQ values, when the 
profile drag and other fundamental characteristics of 
the airfoil section are taken into account, by means of 
simple formulas involving parameters given for each 
airfoil section in a table of airfoil characteristics. This 
method has an important advantage in that the induced 
drag, which is really only a component of the c ?0 at 
each section, is entirely eliminated from the analysis. 

For the purpose of determining the lift distribution 
corresponding to the c iQ values along the span, the lift 
load along the span is considered as being made up of 
two independent parts that will be referred to as the 
"basic lift distribution" and the "additional lift dis- 
tribution." The basic lift distribution is represented 
by the e? 0 distribution along the span when the total 
wing lift is zero. This basic lift distribution, which is 
the distribution arising by virtue of aerodynamic twist, 
may be considered to exist unaltered as the lift and angle 
of attack are changed. The additional lift distribution, 
as the name implies, represents the distribution of addi- 
tional lift associated with changing the angle of attack. 
Wing theory indicates that, as long as the airfoil sec- 
tions of the wing are working within a range of normal 
lift-curve slope, the form of the additional lift distribu- 
tion is the same at all lift coefficients and is independent 
of wing twist, of aileron or flap displacements, and of 
other characteristics that affect only the basic lift dis- 
tribution. Experiment shows that this deduction is 
approximately correct for wings with well-rounded 
tips. For such wings, the additional lift distribution is 
given as a function of the plan form and aspect ratio 
in terms of the additional lift coefficients Ci aV that is, 
the section additional lift coefficients for a wing lift 
coefficient of 1. The lift distribution for any wing is 
then found in terms of the wing lift coefficient G L} the 
basic lift coefficient c lb , and. the additional lift coeffi- 
cient Ci al 

ci 0 =Ci b +C L c lal (1) 
GENERAL PROCEDURE 

MONOPLANE 

It is advisable first to choose a backward fore-and-aft 
reference axis x usually parallel to the reference axis, or 
thrust line, of the airplane and an upward z axis per- 
pendicular to it. (See fig. 1.) Upward and backward 
air forces and distances are thus considered positive. 
Air-force components along these axes are then ex- 
pressed at each section of the wing by 

dX=c x g c dy (2) 



and dZ=c z q c dy (3) 

where X and Z are the components of air load along the 
axes, and c x and c z are determined from c to and the 
known characteristics and attitude of each airfoil sec- 
tion. The pitching moment about the origin contrib- 
uted by each section is 

dM=c maCm q c 2 dy^c x gc z dy—c z qcxdy (4) 

where x and z are distances measured from the origin 
to the aerodynamic center of the airfoil section (see 
table I and appendix B) and the signs of the terms are 
so taken that stalling moments are positive. 




Figure 1.— Airplane drawing and balance diagram. 

Thus far the origin has been arbitrarily chosen. If, 
with this arbitrarily chosen origin, the coordinates x a . c . 
and z a .c. of the aerodynamic center of the entire wing 
(fig. 1) are found, the origin of coordinates may then 
be moved to this point and from equation (4) there 
may be determined a value of M a .c.!q_ that has sensibly 
the same value for all flight conditions. 

Aerodynamic center and additional lift distribution. — 
For the purpose of finding the aerodynamic center of 
the wing, it is necessary to consider only the additional 
distribution. In fact, the aerodynamic center of the 
wing may be considered as the centroid of all the addi- 
tional loads. For wings with linear taper and rounded 
tips, values of L a , giving the load distribution for C L =1, 
may be found from table II for various sections along 
the span. The values of L a were derived as outlined in 
reference 2. The corresponding values of Ci al for various 
sections along the span are found from the relation 

c 'ai = ~X* ^he c o rres P onc Jing values of c dQ at each 
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section are calculated using the method indicated in 
figure 2 or, if the profile-drag polar curve for the section 
is available, they may be read from it. Then 

c *ai = c <i 0 cos 6 Za - c , al sin B Za (5) 

and 

C *al = C lul COS d !a + C « 0 Sil1 K ( 6 ) 



in which 0, 



'hi 



4 oti — i; do, the section lift-curve 



slope, and a to , the angle of attack of zero lift, are given 
in table I; and i is the incidence of the chord at each 
section with respect to the x axis. 



.036 
.032 
.028 
1.024 
' .020 
" .0/6 



.004 



0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 

\ C l~ C l opl\/ C l ,.a " C lopl 

Fxguke 2.— Chart for the determination of the profile-drag coefficient cj 0 at 
any lift coefficient ct 0 . 

The next step is to plot c Xal c, c Zal c, c X(ll cz, and — c Z(ll cx 
against y and to fair curves through the plotted points. 
Twice the area under each curve from y=Q to y~b/2 
is then, respectively: X al /#, Z al /g, M Xal /q, M z Jq. The 
coordinates of the aerodynamic center of the wing are 
then found 

—M Zgl 




M 



x„ 



(8) 



Pitching moment about the wing aerodynamic cen- 
ter. — The additional load distribution for 0 L =l and 
the position of the aerodynamic center are now known. 
The next step is the determination of the basic load dis- 
tribution (that corresponding to C L =0) and from it the 
basic pitching moment or the aerodynamic pitching 
moment of the wing about the aerodynamic center. 
The basic distribution for wings with linear twist may 
be obtained from table III in terms of the load pa- 
rameter L b for a number of sections along the span. 
The method of deriving the L b values is given in refer- 
ence 2. When the wing has partial-span flaps, the 
basic distribution may be obtained from reference 3, 



Following the system that was previously used, c lb 
values corresponding to the basic lift distribution are 

found for each section from c h =L b ^-, corresponding 
Ca 0 values determined, and c Xb and c Zb calculated from the 
formulas 

c n ^c do cos d H — Ci b sin 6 ib (9) 



Vz b = c h cos 6 Zb +c do sin 6 Z 



(10) 



where 0 Zft =— + a to —i. 
a 0 

Likewise are plotted curves of c Xb cz' and ~c Zb cx', 
where ?J and x' are the now coordinates of the section 
aerodynamic center from the aerodynamic center of 
the wing. The areas are then determined. In addi- 
tion, another curve formed by plotting c m(i c c 2 is drawn 
and the area determined. Twice these areas then give, 
respectively, (M x J<i) ac} (.M z Jq) aC) and M s /q. The 
desired wing pitching moment about the aerodynamic 
center is found from 



(11) 



lift distribution and total lift. — When the total wing 
lift or normal-force coefficients are known or specified 
by design conditions, the force distribution may be 
found immediately in terms of ci Q values along the 
span from 

C 'lQ — c '6+^ c f«i 

For wings having well-rounded tips, the lift distribution 
may thus be found in terms of the c h and Ci at values 
previously determined. This method will give a good 
approximation of the actual lift distribution in such 
cases. When, for any reason, the tip loads are of 
critical importance, that is, if the wing is tapered less 
than 2:1 and has a tip blunter than semicircular, the 
lift distribution should be determined according to the 
method given in appendix A or reference 4. If the 
wing plan form departs from a straight taper, the lift 
distribution should be determined from suitable theo- 
retical methods (references 2 and 3). In any event, 
the loads are represented by the Ci 0 distribution and 
may then be resolved to give chord and beam compo- 
nents and moments. 

In general, the wing lift coefficient GY for the steady- 
flight condition preceding an accelerated-flight condi- 
tion will be first determined. After the tail load and, 
finally, the wing lift L are determined from the balance 
diagram for the steady-flight condition, the correspond- 
ing wing lift coefficient is found 



(12) 
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The wing lift coefficient C L for an accelerated-flight 
condition may then be determined. Tor example, it 
may be that the acceleration is produced by a sharp- 
edge gust, and the wing lift coefficient is determined 
by the simplified formula 

C L =C L '+m% (13) 

where 6Y has just been found, U/V is the ratio of the 
gust velocity to the flight velocity, and m is the slope 
of the wing lift curve, which may be determined from 
the values of a 0 or m 6 , tabulated for the airfoil sections, 
by employing the method indicated later in figure 12. 

The required lift distribution is then found in terms 
of the value of c io at each section from 

Ci Q =c lb +C L c lal (14) 

From these values of the lift coefficient at each section, 
the required coefficients representing the components 
of the air load may be computed and the total load 
components then determined as before by measuring 
the areas under curves representing c x qc and c z qc. 
Some question exists, however, in regard to the values 
of c do that should be used in the computation of c x and 
c 2 for the accelerated-flight condition. 

Conditions and forces encountered instantaneously in 
accelerated-flight conditions after a suddenly changed 
angle of attack. — In an accelerated-flight condition the 
C L value calculated from (13) and the c }o values from 
(14) may exceed the maximum lift coefficients. Such 
conditions are possible on entering a sharp-edge gust or 
in abrupt maneuvers owing to the considerable time 
required to accumulate the increased volume of re- 
duced-energy air associated with the increased bound- 
ary-layer thickness or separated flow that will finally 
prevail at the increased angle of attack. Lift values 
should be based on the lift-curve slope extended without 
regard to the usual burbling. Such lift values are ob- 
tained simply by following the outlined procedure. 
The c do values, however, deserve special consideration. 
The increasing profile-drag coefficients at the higher 
lift coefficients are likewise associated with a thickening 
boundary layer or a separating flow that will not occur 
at once when the angle of attack is suddenly increased. 

The profile-drag coefficient for these transient condi- 
tions for a given lift, whether or not the lift exceeds the 
value given by wind-tunnel tests as the maximum, is 
undoubtedly less than the profile drag determined in 
the wind tunnel under steady conditions. The true 
value, however, is unknown and, in fact, a series of 
values increasing with time will exist. It may there- 
fore be expedient in some cases to determine the force 
components on the wing by assuming that c do retain its 
initial steady-flight value throughout the subsequent 
relative pitching motion of the wing. On the other 
hand, if it is desired to investigate the higher values that 



the profile drag will later assume, c d(] may be found in 
the usual way from the wind-tunnel data unless c io is 
greater than Ci max , in which case some value of c dQ may 
be assumed. The value c ao —0.l is suggested. 

The distribution of the resolved components and 
moments and the total wing components. — Values of 
c dQ and Ci 0 for the sections along the span having now 
been established, the distribution of the air-force com- 
ponents, given by values of c x and c z , may be found from 

c x =c<i 0 cos 6 Z — c, 0 sin 6 Z (15) 
c z = c, 0 cos flz + c^ sin d z (16) 

where 

The torsional moment contributed by each section 
about its aerodynamic center is simply 

dM Sa c =c ma c qc 2 dy (17) 

For some problems, components and moments with 
respect to axes in the wing may be desired rather than 
the components given by c x and c z with respect to the 
airplane. For example, "chord-truss" and "beam" 
components may be desired at each section. These 
components represented by c c and. c b may be obtained 
from a slight modification of (15) and (16). 

c c — Cd 0 (l +tan d c tan 4>) cos 6 C (18) 
— Ct 0 (l— cot 6 C tan 4>) sin 0 C 

c b =Ci 0 (l — tan 6 b tan <£) cos d b (19) 
+ c rfo (l+cot Ob fan cj>) sin d b 

where 

a Cl " i 

a C ' h I 

Ob^---TOCi 0 — l b 

(j)=i b —i c 

and i c is the incidence of the section chord with respect 
to the chord-truss direction (plane of the drag truss) 
and i b is the incidence of the section chord with respect 
to the perpendicular to the beam direction (the per- 
pendicular to the spar web). The distribution of the 
chord and beam components C and B may then be 
calculated from 

dC—c c . qc dy (20) 

dB=c b q c dy (21) 

Torsional moments contributed by the sections about 
some axes in the wing other than the axes of the aero- 
dynamic centers of the sections as, for example, the 
wing torsional axis, may be desired in some instances. 
The moment about the torsional axis M T is found from 

dM T =c mac q c 2 dy+c c qc z T dy+c b q c x T dy (22) 
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where z r is the distance of the torsional axis below the 
chord plane through the aerodynamic center and x T 
is the distance of the torsional axis behind the beam 
plane through the aerodynamic center of the airfoil 
section. 

The total forces and moments may then be found 
from the components or, more conveniently for Z and 
M a .c, from, the summations previously made: 



Z Z^ . n 



z_ a 



and Ma.e./d is a value obtained from (11). In order to 
find X/q, however, the c x components should be summed. 

Permissible approximations. — -For the practical ap- 
plication of this method, certain approximations will 
often be justifiable. The approximations that will be 
found convenient and usually justifiable are made by 
assuming that 

cos 02=1 

and 

c<i 0 sin 03 = 0 

The magnitude, but not the direction, of C L and G N 
may then be taken as the same; the following quantities 
are also equal in magnitude but not in direction: 

BIPLANE 

The present unsatisfactory status of biplane theory 
and the large number of variables in the biplane shape 
or arrangement combine to prevent a completely ra- 
tional solution of biplane problems by either theoretical 
or empirical methods. It is possible, however, to 
compute the forces and moments on "conventional" 
biplane wings by scmiempirical methods that give fairly 
satisfactory results. 

In general, the biplane calculations follow the prin- 
ciples and procedure previously outlined for the mono- 
plane, the main extensions therefrom lying in the de- 
termination of the lift distribution between the wings 
and the determination of the biplane effect on the mo- 
ments of the individual wings. The lift distribution 
between the wings is found according to the method 
developed by Diohl in references 5 and 6; the biplane 
effect on the moments of the individual wings is found 
according to a procedure outlined later in this report. 
Although a biplane has no aerodynamic center, a locus 
of points about which the pitching-moment coefficient 
of the cellule remains constant can be found. This 
locus is analogous to the aerodynamic center of the 
monoplane but lacks its practical utility. Neverthe- 
less, since it leads to a better understanding of biplane 
phenomena, the locus of points of constant pitching 
moment will first be discussed. 

Locus of points of constant pitching moment. — Ac- 
cording to DiehPs solution of the lift distribution be- 
tween the wings, the lift coefficients of the individual 



wings plotted as functions of the biplane coefficient are 
straight lines that intersect at some value of the biplane 
lift which is, in general, not equal to zero. A typical 




Figure 3.— Typical biplane wing lift curves. 



Y ax is 




Figure 4.— Force diagram for determination of* and z. 

case is shown in figure 3. These individual wing lifts 
may be considered to have their points of application at 
the aerodynamic centers of the individual wings, be- 
cause, as will be indicated later, the monoplane value of 
the aerodynamic center of either wing is not affected by 
the opposite wing; only the basic moment is affected. 

Now, if it be assumed that the biplane lift relations 
are equally applicable to the Z components, 1 it is clear 
that the location of the center of the Z components may 
be considered fixed in the direction of x, the ratio of the 
change in Z force on the upper wing to the change on the 
lower wing being constant. Reference to figure 4 
shows that the x location of the locus of constant mo- 
ment can be found from the relation 



— JC- 



+Sijx L 



1 This assumption is perfectly valid in this case, since the slight error involved is 
within the error of the scmiempirical method of determining the lift distribution. 
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in which K 2 is Diehl's biplane lift function, as indicated 
in figure 3. 

Unlike the ratio of the Z forces, the ratio of the X 
components is not independent of the biplane lift be- 
cause of the nonlinear relation between profile drag and 
lift in combination with the inequality in lift on the 
upper and lower wings and because of the trigonometric 
relation between the lift and its X component. The 
point about which the pitching moment remains con- 
stant therefore moves in the z direction with changes in 
lift or in X force ratio. Thus, according to figure 4, at 
any value of the biplane lift for which the X components 
may be determined 

- X L Q 

Xjj-\rX L 

A graphic illustration of the behavior of z is given 
in figure 5, which shows values calculated for the bi- 
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Figure 5.— Variation of z wUSi biplane lift coefficient. 

plane selected for the illustrative example given later 
in the report. At the higher values of the lift coefficient, 
the points of constant moment are close to the upper 
wing. In this condition both upper and lower X 
components act forward, the upper component being 
the larger. At a lift coefficient of about 0.33, z is 
indeterminate because the upper and lower X compo- 
nents are equal in magnitude but opposite in direction. 
In this condition the resultant force is in the z direction 
and the X components form a pure couple. At the 
lower lift coefficients both X components act rearward 
and are of nearly equal magnitude so that z is approxi- 
mately half the gap. 

It can be seen from the foregoing discussion that the 
biplane has no useful counterpart of the monoplane 
aerodynamic center. For this reason, biplane problems 
are best solved by proceeding directly to a solution of 
the forces and moments. 

Lift coefficients of individual wings. — The first step 
in the biplane solution is to determine the lift coefficients 
of the individual wings as functions of the lift coefficient 
of the cellule. As previously indicated, this step may 



be performed according to the method developed by 
Diehl in references 5 and 6. When this method is 
used, however, it is recommended that, in cases involv- 
ing large negative stagger, values of K 2 q be determined 
from a curve faired through the experimental points of 
figure 13 of reference 5, rather than from the linear 
relation between K 20 and s/c (equation (15a), reference 
5). 

Distribution of force components.— The wing lifts 
corresponding to any biplane lift having been found, the 
force distribution on the individual wings is determined 
in the same manner as for monoplanes. This procedure 
neglects the effect of interaction of the individual 
wings and leads to some error, which is probably small 
in practical cases. 

Pitching moment of biplane cellule. — The pitching 
moment of the whole cellule about any arbitrarily 
selected Y axis is found in the same manner as for the 
monoplane from a summation of the moments due to 
the Z and X components of force and to the section 
characteristics. To this total moment a correction, 
constant throughout the lift range, may be applied to 
staggered arrangements to obtain a more accurate 
result. 

The correction is based on the fact, indicated by 
available test data, that the couple created by the lift 
forces on the individual wings of a staggered biplane 
with no decalagc at zero cellule lift is exactly balanced by 
predominating increments of moment on the individual 
wings plus a secondary couple due to the biplane effect 
on the drags of the individual wings. The moment 
correction, therefore, constitutes simply a subtraction 
of the couple created by the Ki forces due to thickness- 
gap ratio, stagger, and overhang from the total moment 
M previously found. Thus 

M Yr = 7iM-(K w +Kn+K u ) Su sg_ 

where K i() , K n , and K n are Diehl's lift functions for 
thickness-gap ratio, stagger, and overhang and s is the 
stagger measured between the aerodynamic centers of 
the individual wings. 

The function K 12 , which is due to decalage, is not 
included in the correction. 

Pitching moments of individual wings. — As previously 
mentioned, the couple due to the Ki forces, if decalage 
is neglected, is exactly balanced by predominating op- 
posite moments on the individual wings and a less im- 
portant couple due to biplane effect on the drags. This 
drag moment is small compared with the Ki couple and 
therefore negligible, since the K x couple itself is small. 
The Ki couple may therefore be considered to be en- 
tirely balanced by increments of moment on the indi- 
vidual wings. No information exists, however, as to 
the distribution of these moment increments between 
the upper and lower wings ; a consideration of this prob- 
lem led to the conclusion that a reasonable assumption 
would be to divide the balancing couple equally between 
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the wings. This assumption leads to very low incre- 
ments of pitching-moment coefficient on the individual 
wings; in several cases that have been examined the 
values were well below 0.01. In view of such low values 
and the uncertainties in regard to the distribution, it is 
believed advisable to neglect these increments in com- 
puting the pitching moments of the individual wings. 

Another biplane effect on the individual wing mo- 
ments, however, should be taken into account. Its 
physical cause is not known at present, but it is probably 
due to the profile drag of the wings, which results in a 
pressure gradient from the leading to the trailing edge 
between the wings and to the curvature in the stream- 
lines at each wing induced by the opposite wing. An 
examination of test data obtained both in flight and in 
wind tunnels showed that this biplane effect on the 
wing moments is, for all practical purposes, a linear 
function of the thickness-gap ratio given by the relation 

Acmo^ =0.l£ 

These increments, for the data available, do not notice- 
ably contribute to the resultant biplane moment; the 
total increment of moment on the upper wing must 
therefore be approximately equal and opposite to that 
on the lower wing. 

In order to effect the practical application of these 
increments to the wings, it is assumed: (1) That the 
increment is distributed along the entire span of the 
shorter wing but only along that portion of the span of 
the longer wing that lies within the projected span of 
the shorter wing; and (2) that the increment of pitching- 
moment coefficient is distributed uniformly along the 
span of each wing between the limits of the pitching- 
moment distribution. On the basis of assumption (1), 
the value of hc mo . t . is found for the upper wing from 

the foregoing relation using the average value of t/G 
based on the lower wing for the portion of the span 
affected. Then 

Ar — — Ar \y^u\yCu' 

where S v ' is the area of the portion of the upper wing- 
involved. 

S L ' ', the area of the portion of the lower wing in- 
volved. 

c v f , average chord of the portion of the upper 

wing involved. 
c L ', average chord of the portion of the lower 

wing involved. 

LOAD DISTRIBUTION OVER AIRFOIL SECTION 

The solution of the general problem has been com- 
pleted except that the distribution of the air forces 
along the chord at each section has not been determined, 



the net section lift, drag, and pitching-moment coeffi- 
cients having been employed heretofore rather than the 
distributed air loads at each section. Although the 
distribution of the air load around the airfoil section 
may not always be required, this distribution will be 
considered in order to make the analysis complete. 

General procedure. — The previous analysis gives the 
section lift coefficient Ci 0 , the method of finding the 
normal- and chord-force coefficients c n and c c , and the 
pitching-moment coefficient c ma c _ at each section cor- 
responding to any given loading condition of the com- 
plete airplane with which the designer is concerned. 
The corresponding distribution of the air load over the 
section will be given in terms of the normal-force co- 
efficient by giving the distribution of the normal-pres- 
sure coefficient P along the chord of the section. Of 
course, this distribution gives no chord force but the 
chord force is known and may be considered as applied 
at the aerodynamic center. Its distribution will not 
be considered, the chord force being small and dis- 
tributed over only a small distance equal to the wing 
thickness. Although the moment contributed by this 
distribution cannot be entirely neglected, the normal- 
force distribution will be slightly modified, more or less 
arbitrarily, so that it will give exactly the correct pitch- 
ing moment about the aerodynamic center. 

Determination of normal-pressure coefficients. — As 
previously stated, the distribution of the air load along 
the chord is found by determining the normal-pressure 
coefficient P, that is, the ratio of the pressure difference 
that may be considered as acting at any point along the 
chord to the dynamic pressure q. The distribution is 
defined by the values of P at a number of points along 
the chord. As with the span load distributions, it is 
convenient to consider the distribution as made up of 
two independent parts, one the distribution for zero 
normal force P 0 and the other an additional distribu- 
tion giving all the normal force. The total normal- 
pressure coefficient at each point is then 

■ P=P 0 +c n P a (23) 
The value of P a is found from 

P a =P al +^P a . c . (24) 

where values of P aX and P a . c . are given by curves and 
tables for typical airfoils in figure 6. The designation 
of the airfoil class in this respect corresponds to a letter 
given for each section in the PD column of table I. 
Values of x a . c ./c are also found from table I by dividing 
by 100 the x coordinate of the aerodynamic center. 
A single table of the P a . c . distribution, which is taken 
as the same for all airfoils, is given in figure 6 . 

The value of P 0 is found from the so-called "basic 
distribution," thus 

Po=P*-c nb P a (25) 
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Note.— Type A distributions have not yet been determined. 

Figure 6. — Pressure distribution— additional. 

The basic distribution P h and the basic normal-force 
coefficient c nb are, in turn, found as the sum of two parts 
due respectively to moment and camber, thus 

P*=-c ma . c P> m +jP b c (26) 

and 

Cny=— Cm a . c .(Cn) bm -\-f(fi n ) be (27) 

Values of P bm and the corresponding values of (c n ) bm are 
given in figure 7, as well as values of P hc and (c n ) bc for 
airfoils of classes as indicated in the airfoil table by the 
number following the letter in the PD column. For 
example, the number 10 indicates that P im is class 1 
and P bc is class 0. The zero signifies that P bc and (c n ) u 
are both zero. The values of c ma c and the section 
camber z c fc are both found from table I, zjc being 
found by dividing the mean camber as given, in percent 
of the chord, by 100. 



When the actual calculation for any given airfoil 
section is made, values of P 0 and P a should be calculated 
and tabulated for the standard stations along the chord. 
For any section normal-force coefficient c n , the corres- 
ponding values of P are then found simply from (23) 
by multiplying the values of P a by c n and adding to P 0 - 
The actual pressure difference acting at each point in 
pounds per square foot is, of course, obtained by multi- 
plying by the dynamic pressure in consistent units. 
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Figure 7. — Pressure distribution— basic. 

Finally, consider briefly how the air pressures are 
divided between the upper and the lower surfaces. 
Pressure-distribution diagrams, given elsewhere, indi- 
cate the pressure on the upper and on the lower surface 
as measured from the static pressure as a reference. 
The designer, however, is not primarily concerned with 
these pressures but with the pressure differences across 
the wing covering which, of course, produce the air load 
on it. These pressure differences are a function of the 
internal pressure, that is, the pressure within the wing. 
If the wing is well vented, the internal pressure and the 
upper and lower covering loads may be estimated. For 
this purpose the lower-surface pressure distribution is 
estimated, remembering that the positive pressure can- 
not exceed by lg the static pressure, and the upper- 
surface distribution determined from the known values 
of the differential-pressure coefficient P. If greater 
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accuracy is required, the method of reference 7 or the 
results of reference 8 may be employed to calculate the 
pressure distribution on the lower surface. 

SAMPLE CALCULATION 
MONOPLANE 

In order to make this example as general as possible, 
a case is chosen for which the design condition repre- 
senting a 30-foot-per-second gust encountered at high 
speed causes the lift coefficient to exceed the usual 
maximum lift coefficient for the airfoil. The example 
does not, however, deal specifically with the procedure 
to be followed in cases for which portions of the wing 
are replaced by the fuselage or nacelles. The treat- 
ment, nevertheless, is exactly the same in such cases if 
the standard N. A. C. A. wing area, including those por- 
tions of the wing imagined as inside the fuselage or 
nacelles, is used for S. The solution is thus found by 
considering those portions of the wing to be actually 
present and undisturbed, the wing being imagined as 
extending continuously over those portions of the span. 
The calculated loads for these imaginary portions of the 
wings may later be applied to the fuselage and nacelles. 
In extreme cases a special treatment may be required. 
A wing of the U. S. A. 35 type is chosen so that some 
aerodynamic twist will be present in spite of the fact 
that the wing is not twisted with reference to the air- 
foil chords. The drag truss, for generality, is taken at 
an angle to the plane of the airfoil chords. The 
analysis is begun from the airplane drawing in figure 1 
and from the following data on the airplane and wing: 



Weight 1,0001b. 

Power 35 hp. 

Propeller efficiency 75 percent. 

High speed 95.3 f. p. s. 

Wing incidence 4°. 



Wing: U. S. A. 35 type, aspect ratio 5, rounded tips, 
area 180 sq. ft., root chord 8.268 ft., taper ratio 
0.5, no geometric twist, beam direction perpen- 
dicular to chord, drag truss (chord direction) in- 
clined upward at the leading edge with respect to 
the chord 4° at root to 2° at tip, airfoil section at 
root U. S. A. 35-A, at tip U. S. A. 35-B. 

Calculation of wing aerodynamic center—The first 
step in the procedure is to choose the reference axes. 
The axes are chosen, for generality, originating at the 
center of gravity with one axis parallel to the thrust line 
although, in this instance, some simplification might 
have resulted from choosing an axis in the direction of 
an airfoil chord because this direction is the same 
along the wing (no geometric twist) and perpendicular 
to the beam direction. Table IV is then filled in to 
give the necessary data for computing the aerodynamic 
center of the wing. The various columns leading first 
to the calculation of the additional-load curves for 
G L ~1 and finally to the position of the wing aerody- 
namic center are filled in as follows : 



Column 1— Stations along the span chosen arbi- 
trarily or to agree with those in table II. These sta- 
tions are indicated on the airplane drawing (fig. 1). 

Column 2. — Values of L a from additional-load table 
(table II) for aspect ratio 5, taper ratio 0.5. 

Column 8. — Values of c from the airplane drawing. 

Column 4- — Values of c ial from the multiplication of 
(2)bySJcb. 

Column 5. — Values of a 0 from airfoil characteristics 
(table I) interpolating between U. S. A. 35-A and U. S. A. 
35-B sections for intermediate sections of wing. 

Column -Values of c lal /a 0 from (4) and (5). 

Column 7. — Values of oti 0 by the same method as (5). 

Column 8 — Values of — i, the incidence of the wing 
chords with respect to the x axis with reversed sign, from 
airplane drawing. 

Column 9. — Values of c do . The profile-drag coeffi- 
cients are calculated for each section as indicated in 
table IV-A. The thickness ratio of each section t/c is 
obtained from the airplane drawing. Minimum pro- 
file-drag coefficients c da are obtained from a curve of 

min 

profile-drag coefficient against section thickness, paral- 
leling the typical curve given in reference 9 (fig. 91) but 
passing through the values indicated in table I for the 
U. S. A. 35-A and U. S. A. 35-B sections. Values of Ci opt 
and Ci max are obtained from table I. From the values 

in the preceding columns, the ratio \ Clal . — is com- 

puted. From this ratio and the curve of figure 2, the 
Ac d(s values are obtained, which are added to the values 
of c d( . to give the desired c dn . 

min 

Column 10. — Values of 6 Za . From the addition of 

(6), (7), and (8), where 9, a is (^-fa^-i). (See equa- 

\ a 0 / 

tions (5) and (6).) 

Columns 11 to 16. — From preceding columns. 

Column 17. — -Values of c Xal from (13) + (14) following 
equation (5). 

Column 18. — Values of c, al from (15) + (16) following 
equation (6). 

Column 19. — Values of z from the airplane drawing, 
upward coordinate of aerodynamic center of section. 
May be obtained from airplane drawing after locating 
the aerodynamic center of the tip and the center sec- 
tions from table I. The corresponding aerodynamic- 
center positions for the intermediate sections may be 
taken along the straight line joining these points except 
for the rounded-tip sections. 

Column 20. — Values of x, backward coordinate of 
aerodynamic center of section, obtained from the air- 
plane drawing as with (19). 

Columns 21 to 24- — Products from previous columns. 

These pitching-moment and loading results are 
plotted as in figures 8 and 9, and the areas measured to 
find M Zal !q, M Xal /q., Z a i/g_, and X a i/q. From these values 
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the coordinates of the wing aerodynamic center are 
found from equations (7) and (8). 
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Figure 8.— Plots for the determination of the components of the additional 
wing pitching moments and the wing aerodynamic center. 
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Figure 9 — Distribution of additional z and x components and determination 
of additional wing components. 

Calculation of wing pitching moment about aerody- 
namic center. — The next step is to carry out practically 
the same procedure for the basic load distribution in 
order to find the wing pitching moment about the aero- 
dynamic center. The origin of coordinates is moved to 



the aerodynamic center and another set of calculations 
is made for the basic load distribution as indicated in 
table V. The only differences worth noting are the 
different values taken from the tables, values for the 
basic load distribution L b from table III, and the method 
of obtaining from these the c lb values in table V. The 
Ci h values follow from those in the second column taken 
from table III, by multiplying by ea^Sjcb. 

The term ea 0 takes into account the aerodynamic 
twist of the wing, which is assumed to vary linearly 
along the span. The twist e is measured with respect 
to the zero lift directions for the center and tip sections, 
being positive when the effective incidence is washed in 
from the center toward the tip. It is evident that the 
term ea 0 is a Ci difference between airfoil sections corre- 
sponding to the center and tip sections when the section 
angles of attack have the same relation as in the wing. 
In other words, ea 0 may be calculated as follows : 

ea 0 = [d 0 (ol — a , 0 ) ] t i P — [a Q {a—Oi iQ ) ] cen tar 

This procedure is strictly correct theoretically only 
when a 0 does not vary along the span. When a 0 
varies, the best practical result is probably obtained by 
calculating ea 0 as a Aci for an a near the value at which 
the load distribution is desired. 

The value of a may be taken as zero for the center 
section and, because no geometrical twist is present, 
the value of a at the tip is then also zero in this instance. 



edo — ( — &oO:; 0 ) tip ( ttoaz 0 )cen/er 

= [- (0.099) (-5.2)]-[- (0.095) (- 
= 0.515-0.760^-0.245 



5.0)1 



Values of the factor ea 0 S/cb are then obtained at 
each station along the span by which the values taken 
from table III are multiplied to obtain the c* 6 values. 
From the c ib values, the calculations proceed to the final 
res^TsT^nich are given in the last columns of table V. 
These results are plotted and the areas determined to 
find (M x Jq) ax ., (M z Jg) a .c., and (M s /q) a . c . as in figure 
10. These values are added to obtain M fliC ./g- 



M a . c 



=0+4.84— 113.36- -108.5 



which, multiplied by a, gives M a .c, the required pitch- 
ing moment of the wing about its aerodynamic center. 

Calculation of forces and moments in accelerated- 
flight condition.— The exact procedure to be followed 
from this point on is dependent on the result desired. 
If a result meeting arbitrary design requirements is 
desired, the particular specified procedure will be fol- 
lowed. If, on the other hand, the most reliable actual 
air loads for a given design condition are desired, 
another procedure may be advisable. 

From the method of references 1 and 10, for example, 
the applied load factor Ui is determined and the wing 
normal-force coefficient C Nl is taken as n^/q^ where 
s is the effective wing loading and q L is the dynamic 
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pressure for the design speed. Corresponding values 
of the chord-force coefficient c c are obtained as more or 
less arbitrarily specified, and the pitching character- 
istics of the wing are rather arbitrarily given by speci- 
fying that the center-of-pressure position be taken as 
the most forward position for the wing between 
G L =-G Nl and C L =C Lmax) unless C Nl exceeds C Lmax , in 
which case a value taken from the extended center-of- 
pressure curve is specified. After G L is calculated from 
the specified C Nv the lift-coefficient distribution may be 
found by adding the basic and additional lift coeffi- 
cients in accordance with the relation 

and including, when necessary, the tip corrections 
given in appendix A. The corresponding specified 
values of the center of pressure and of c c may then be 
applied at each section and the forces and moments 
resolved as desired* for structural analysis. 



y.tt. 



-2 



-4 



-6 



-7 











1 

h cx 


' Ar 


ea - 


= 2. 


42 f 


t- 3 (j 




184 


i 

ft. 3 


























— X 


























cz \ 


























Area = 


0 
































y 
































































































































































































c 


m a.C. 




























Area = - 
i i 


56. t 


"6 ft 


s 
























36 ft 










































7 

























































































Figure 10.— Plots for the determination of the components of the basic wing pitch- 
ing moment. 

The foregoing procedure, however, will not be fol- 
lowed in this example. Specified design conditions and 
methods vary and, in many instances, it is believed 
that designers will wish to investigate loadings under 
conditions other than those specified. The example 
will therefore be carried through using the procedure 



that may be expected to give the best approximation 
to the actual air forces. 

The first step is to obtain the lift coefficient G L r 
corresponding to the steady-flight condition before 
entry into the gust. For the present example, this 
condition is represented by high-speed level flight. 
The corresponding 6V value is obtained from the 
balance diagram for the airplane for this condition. 

For the construction of the balance diagram, it is 
necessary to know the angle of the flight path so that 
the direction of the weight vector may be determined. 
A trial value of C L f is first taken, assuming a down tail 

load of 30 pounds, — ^ — =0.530. The wing angle of 

attack as measured by a s , the angle of attack referred 
to the chord of the center section, may then be deter- 
mined by the method indicated in figure 11: 
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Determine the angle of attack from: 

<* 3 = Cj,/a+ (<*[,),+ Jt 

where a„ angle of attack referred to the chord of the central section of the wing. 
Cl, wing lift coefficient. 
a, wing lift-curve slope per degree. 
( a 'o).» ar) gl e 01 zero lift 01 tfrfc central section. 

e, angle of aerodynamic twist. 
A, aspect ratio. 

CL=a[a.—(<x] a ) t —Je] or, angle of zero lift for the wing referred to the chord 
of the center section=(ai 0 ) l - r -Je. 

Figure 11.— Lift in terms of angle of attack for tapered wings with twist. 

From figure 12 

a=f 



a 0 



57.3a 0 
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Taking a 0 —0. 
the wing 



as a mean value for the sections of 



a=0.999- 



0.096 



=0.0711 



57.3X0.096 



7f5 



The angle of zero lift for the root section (a jQ ) t is taken 
from table I. The twist e is ea 0 /<z 0 or -0.245/0.096= 
— 2.55°. The factor J from figure 11 is —0.408. 
Then 



a 

0.530 



0.0711 

=0.5° 



■8 + (-0.408) (-2.55) 



As the incidence of the center section is 4°, the angle 
of the thrust line with the horizontal is 0.5°— 4.0°= 
—3.5°. The weight vector may therefore be drawn as 
indicated in figure 1 and the pitching moments may be 
taken about the wing aerodynamic center to determine 
the tail load F t . The dynamic pressure is 

g=I (0.002378) (95.3) 2 
= 10.79 lb./sq. ft. 

= -108.52X10.79 
= -1,171 lb.-ft. 

Although for other purposes, such as balance calcula- 
tions, a better moment analysis may be necessary, the 
thrust moment in this case may be determined with 
sufficient accuracy on the assumption that three-quar- 
ters of the thrust is used in overcoming parasite drag, 
which may be assumed to act approximately along the 
thrust axis and therefore to contribute no moment 
about the wing aerodynamic center. The thrust is 



35X550X0-75 
95.3 



151.5 lb. 



Then writing the moment equation 

(0.25X151.5X3.82) + (1,000X0.82)-1,171-F,14.62=0 
2^= -14.3 lb. 

The final value of the lift coefficient for steady flight 
6Y may then be computed 

1,000+14.3 
Cl 10.79X180 

=0.522 

The new wing lift coefficient C L after entry into the 
gust is now determined from the slope m of the wing 
lift curve and the gust velocity U 



^=0.522+4.07^ = 1.803 

^-Elliptical wing 




6 8 10 12 14 16 18 20 
Aspect ratio, A 

Compute the lift-curve slope from the equation: « 

S7.3«o , Co , j „„„„„-, 

m =f 1+si.wa or a=/ nwJS&A (per degree) 

where m, lift-curve slope (per radian). 

ao, section lift-curve slope (per degree). 
a, wing lift-curve slope (per degree). 
A, aspect ratio (b 2 /S). 
f, plan-form factor. 

When the lift-curve slope is normal, the following approximate equations may be 
used: 

67.3ao , tth 

1+1.843/A 0r m ~ } 0.761+1.403M 

where mt, lift-curve slope (per radian) for wing of aspect ratio 6 with rounded tips, 
laper ratio- c> = c?Ilter chord 
Figure 12.— Values of/ for computing the lift-curve slope. 

where the value of m is found from figure 12 and table I. 

The lift distribution is now determined by calcu- 
lating the Ci Q values in table VI from 

Ci 0 =c h +O L c lal 

where c tb is taken from table V and c laU from table VI. 
The calculations indicated in table VI proceed then to 
the determination of chord and beam components from 
equations (18) and (19). 

It will be noted that the steady-flight value of the 
profile-drag coefficient c da ' has been used in the accel- 
erated-flight condition. This procedure should be fol- 
lowed when a large forward-acting chord force is critical 
for the structure. 

The last three columns of table VI give the required 
data on the air-force distribution as chord and beam 
forces and pitching moments per running foot of span. 
In order to complete the balance diagram, however, the 
total air forces and moments on the wing are required. 
The pitching moment of the entire wing in this case is 
the same as that previously found for the steady-flight 
condition, because M a . c .lg. has not changed. 



= -108.5X10.79=-l,171 lb.-ft. 
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The total wing air-force component Z perpendicular 
to the thrust line may be found with sufficient accuracy 
from the Z a i/<Z value previously determined without the 
necessity of resolving and summing the section forces 

= 1.803X180.0X10.79 
=3,510 lb. 

The total component X for the wing, however, should 
be found by resolving the section forces along the x 
direction and summing to find the total force compo- 
nent X The values of c dQ and c tQ are taken from table 
VI and resolved by equation (15) to find the c x values 
at the various sections along the span. In this example 
a large forward-acting chord component is assumed to 
be conservative so that the profile-drag coefficient 
for the accelerated-flight condition is taken as equal to 
that in the preceding steady-flight condition for the 
determination of c x . These values are then multiplied 
by c, plotted as in figure 13, and the area determined 
to give X. The result is 



X= -42.2X10.79 = -455 lb. 

y, ff - 

O 2 4 G 8 10 12 14 
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Figure 13.— Distribution of X component in the accelerated-flight condition. 



The total air forces and moments are now known and 
may be applied at the aerodynamic center of the wing 
so that the balance diagram may be completed for the 
accelerated-flight condition, thus completing the solu- 
tion with the exception of the determination of the air- 
force distributions over the ribs. 

Calculation of the air-force distribution over a rib. — 
In order to complete the example, the rib distribution 
will be determined for the central section of the wing. 
Reference to table I will show the pressure-distribution 
classification of the U. S. A. 35-A section to be E10. 
The additional pressure distribution is therefore found 
from the class E P ai distribution. Values of P al and 
P 0 . c . are taken from figure 6, and the additional pres- 
sure distribution is then calculated from equation (24) 

p p I x a.e. p 

1 a — * alT c 1 a.c. 

where x a . c . is the distance of the section aerodynamic 
center forward of the quarter-chord point, from table 
I. The calculation may be carried out in tabular form 
as shown in table VII. 

The basic pressure distribution as given by values 
of P b is then found. From equation (26) 

Pb— Cm ac _Pbm~\~~PbC 

The designation in table I of the basic pressure dis- 
tribution for this airfoil section is indicated by the 
number 10 following the E. The designation 10 indi- 
cates that P bm is class 1 and P bc is class zero, that Is, 
P be =0. The basic pressures may then be computed 
as indicated in table VIII, taking values of P bm from 
figure 7 and the value of the pitching-moment coefficient 
c ma e from table I. 

The zero lift distribution given by values of P 0 is 
then obtained by deducting a part of the P a distribu- 
tion corresponding to c Ub according to equation (25) 

Pq — P b Cn b P a 

value of c nb is obtained from equation (27) 

C "b~ ~ C m a .c.( C n) bm -$-~ ( C «) 

v be 

•e the values of (c n ) bm and (c n ) be are given in figure 7, 
he various distributions. In this case (c n ) bm —6.30 
(c n ) b =0, hence 

c W6 =0.111X6.30+0 
=0.699 

F 0 =F & -0.699P a 

or the accelerated-flight condition, the pressure dis- 
ution as given by values of P is then found from 

P = Pq -f" C n P a 

re c„ in this case is the same as the beam-component 



It will be noted that the preceding calculations of 
the forces in the accelerated-flight condition have been 
made on the basis that the dynamic pressure q remains 
unaltered after encountering the gust. This condition 
is possible when the gust has a small angle with the 
vertical, and the procedure is further justified by the 
fact that the gust velocities specified have been largely 
determined on the basis of their effects on airplanes as 
indicated by the simple gust formula without taking 
into account such changes of velocity. In some in- 
stances, however, it may be desired to take into account 
the dynamic-pressure increase due to a gust, in which 
case the gust velocity should not be taken as nearly 
vertical but may be taken at an angle with the hori- 
zontal and the angle determined to give the maximum 
load. 



The 
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coefficient c b and is taken from table VI for the center 
section, 

P=P 0 -fl.695P a 

Finally, the actual pressure differences p are obtained 
by multiplying by the dynamic pressure, g= 10.79 
pounds per square foot. These pressures, calculated 
as indicated in table VIII and giving the final result, 
are plotted in figure 14. 




■20 40 BO 80 100 

Distonce from L.E., percent c 

Figure 14.— Pressure distribution over center section in the accelerated-flight 

condition. 

BIPLANE 

The following example is given to illustrate the appli- 
cation of the method to biplane problems and also to 
illustrate the alternative method of finding the force 
distribution in a case where the empirical tip correc- 
tions may be important. A simple biplane (fig. 15) 









L 4 J - J 




< j 30' — > 



Fuselage ! 

- 



4 '6" 



2'6'~ 



—/?*?' — 




Figure 15.— Biplane cellule for illustrative example. 

has been chosen in order to avoid, as far as possible, 
steps that have already been illustrated in the example 
for the general monoplane. The calculations are made 
for an airplane having the following characteristics: 
Weight 1,636 1b. 

High speed 100 m. p. h. (2=25.6 Ib./sq. ft.). 

Wing cellule: 

Upper: N. A. C. A. 2412 section, span 30 ft., chord 4K ft., 
area 135 sq. ft., no taper, no twist, no dihedral, incidence 
0°. 

Lower: N. A. C. A. 2412 section, span 22 ft., chord 4 ft., 
area (including projection through fuselage) 88 sq. ft. 



Distance from leading edge of upper wing to c. p. of tail 

plane (I), 15 ft. 
Distance from leading edge of upper wing to c. g. (x e . a ), 

1.75 ft. 

The objective of the calculations in this example is 
the solution of the force and moment distribution along 
the upper wing in the 30-foot-per-second gust, corre- 
sponding to design Condition I of reference 1. In 
order to make the example more illustrative, the wing 
lift coefficient for the initial condition of steady flight 
is found from a balance computation, as in the mono- 
plane example, but only the tail load is considered as an 
extraneous force. Because the biplane has no aero- 
dynamic center, an exact balance can be obtained only 
through a process of trial and error; in the example the 
calculations are not repeated to obtain the exact solu- 
tion. 

lift coefficients of individual wings. — The following 
data are pertinent to the solution of the lift distribution 
between the wings: 

Effective stagger s 0 between one-third-chord points 
at zero lift, 1.67 ft. 



Overhang, 



30-22 
30 



=0.267. 



1.67 



= 0.42. 



^=1.125. 



c a 4.5 



t 0.12 



Q 1.125 



=0.1065. 



=0.12. 



With the foregoing data, the method of reference 6 
yields — 

K xo = -0.0178 
i£ n = 0.0123 
K 12 = -0.0274 
K K = -0.0178 

Ki =Xio+Xii+Jfu+Xa = -0.0507 
iPVKao^O.0951 
K n = 0.0083 
i? 22 =:0.0650 

K 2 =F 2 K 20 -\-K 2i +K 22 =0.m4: 
0 Lu = G Lb + {K x +R. 2 C Lb ) 
= 1.168<7£ B -0.0507 



(28) 



C Ll =C Lb — (K 1 -\-K 2 C Lb ) 
=0.742C iB +0.0778 



(29) 



Wing lift coefficient in steady flight (first trial) .- 

Neglecting tail load, 



C T >= 



W 



1,636 _ 



L * (Su+S L )(i (135+88)25.6' 



0.286 
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With this value of C Lb the general method can be 
applied to each wing for the lift coefficients derived 
from equations (28) and (29) to find the moment about 
some arbitrary Faxis and, from this result, the tail load. 
The simplicity of the biplane cellule chosen permits, 
however, a relatively simple solution of the moment. 
Since the wings are rectangular and of constant section, 
the aerodynamic centers of each wing lie on the locus of 
the aerodynamic centers of the sections; the resultant 
wing forces may therefore be considered to apply at 



2 z t 



x 



Ft 




Figure 16.— Skeleton diagram of airplane for biplane example. 

these centers. Reference to figure 16 indicates that the 
moment about the axis 0 may be expressed as 

Z v x„±Z L 2 L +M„+M L +X L G 

in which s is now the stagger between the aerodynamic 
centers of the individual wings, and (Ki—K^Susq is 
the moment correction to allow for the increments of 
moment, which are not taken into account on the 
individual wings. 

For the steady-flight condition the resultant forces 
at the aerodynamic centers of the individual wings are 

Zu =Clu q Su 

= (1.1680^— 0.0507) X25.6X 135 

= 1,035 lb. 
Z L '=C LL 'qS L 

= (0.742<7 Lfl +0.0778) X25.6X88 

= 654 lb. 

The forces X v ' and X L ' are found by summation of 
the force components along the span in accordance with 
the general procedure described in the report. For this 
purpose the span distribution of c, 0 ' (or c z r ) has been 
found according to the alternative method given in 
appendix A, neglecting the tip corrections, which at low 
lift coefficients are very small. 

AV=1.78lb. 
X/=6.14 lb. 

M v = O mo SuCuq = - 684 lb .-ft. 
M L =C mo S h c L <l= -397 lb.-ft. 
K r -K n = -0.0234 
With the foregoing data, and from the airplane 



characteristics, F ? =69 lb., acting downward. The 
corrected value of G Lb is 



lb 



; _jL ,636+69 
25.6 (135 + 88) U,t3U 



which is the value taken for the initial steady-flight 
condition prior to entry into the gust. 

Wing lift coefficient in accelerated flight. — The incre- 
ment in lift coefficient due to the gust is determined 
from the slope m of the cellule lift curve and the gust 
velocity U. The slope of the lift curve may be deter- 
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FiGUKic 17.— Correction factors for wing Cl. 



AO 



mined, as in the case of the monoplane, from the 
expression 

57.3a 0 



57.3oo 
irA 



For the biplane, /may be taken as unity and A— 0 , 0 - 

in which Munk's span factor k may be determined from 
reference 11. In the present example, 



m- 



57.3X0.097 



57.3X0.097 
(1-02X30)* 
*" 135+88 



=3.91 



C Lb — C Lb +m -y 



=0.30+3.91^^=1.10 
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From this point, the distribution of forces and mo- 
ments on the upper wing are to be determined by using 
the method for finding the span load distribution given 
in appendix A. The first step is to find the lift co- 
efficient of the upper wing. 
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Figure 18 —Tip corrections. 



1.0 



1. From C Lb and equation (28) 

G Lu =l. 168X1.10-0.0507=1.233 
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Figure 19— Distribution of forces and moments on upper wing of biplane. 

2. In order to reduce C Lu to allow for the tip cor- 
rection, Fi and F 2 are found from figure 17 

^=0.038 

*W.06 

i^Xi^O.040 

The value of C Lu used to enter the charts is (tables IX 
and X) 

(^/^l^-O.O^^l.lOS 

3. The aspect ratio of the upper wing is 

(30)' 



135 



=6.67 



4. The distribution of c tal is found from table IX, 
which gives the distribution for aspect ratio 6. The 
aspect ratio of the upper wing will be considered herein 
as sufficiently close to 6 to require no interpolation. 
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Values of c tal are tabulated in column 2 and values of 
C Lu "Xci al are tabulated in column 3 of table XI. 

5. Since there is no twist, the values of c ib are zero. 

6. The tip increments are determined from figure 
18 and are tabulated in column 4 of table XI. 

7. The tip corrections are added to the values of 
Ci" to obtain the final C; 0 values tabulated in column 
5 of table XI. 

From this point the procedure follows the general 
method of this report to find the values of c z c of column 
14 and c x c of column 17. The values of Ac do given in 
table XI have been computed for the initial steady- 
flight condition in accordance with the principle of the 
delay in the growth of the boundary layer in accelerated 
flight. These values give the distribution required; 
they are plotted in figure 19. 

In order to find the moment distribution, the basic 
section moment coefficients are tabulated in column 18, 
and to these are added the increments due to tjG and 
the tip effect. The value of Ac mo/ 1 . is found from the 



expression 




= 0.1X0.1065-0.011 

These values are applied only along that portion of 
the span of the wing which lies between the projected 
tips of the lower wing, as indicated in column 19. 

The tip -moment increments are found from figure 18 
and are tabulated in column 20. The resultant dis- 
tribution of c mo is given in column 21, and the final 
values of moment are tabulated in column 22 and 
plotted in figure 19. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., March 25, 1938. 



APPENDIX A 



DETERMINATION OP SPAN LOAD DISTRIBUTION FOR 
SPECIAL CASES 

The tables of span load ordinates (tables II and III) 
referred to in the development of the method and used 
in the monoplane example are, in general, suitable for 
the determination of resultant wing forces and for the 
determination of force distribution for structural appli- 
cations except in cases, such as some externally braced 
wings, in which the tip loading has an important 
influence. For such excepted cases empirical tip cor- 
rections should be applied, in accordance with the 
following procedure. Also, in cases in which the plan 
form departs widely from the straight tapered shape or 
in which there are discontinuities in twist such as occur 
with partial-span flaps, the span load distribution 
should be determined from the basic wing theory. 
For such cases, the method discussed in reference 3 is 
recommended. 

The results from reference 4 are to be used. The 
following procedure should be utilized for obtaining 
the span load distribution with special tip corrections: 

1. From the conditions of the problem determine 
18 



C L , based on the actual wing area, for which the dis- 
tribution is to be determined. 

2. Find AC L {F X XF 2 ) from figure 17, interpolating 
when necessary. Subtract AC L from the value of C L 
found from step (1). 

3. Determine the geometric aspect ratio of the actual 
wing. 

4. From table IX find the Ci al distribution and mul- 
tiply by the value from step (2). 

5. Add to the distribution found in step (4) the c lb 
distribution from table X, reduced or increased in 
proportion to the actual twist. 

6. Find the tip corrections (Acj 0 ) from figure 18(a). 
The affected distance is 40 percent of Sfb. The tip 
increments of figure 18(a) are multiplied by both the 
aspect-ratio and taper-ratio factors given in figure 
18(c). 

7. Add the Aci 0 increments, corrected for aspect 
ratio and taper, to the distribution of step (5). 

8. Add to the section c mo values the tip Ac m incre- 
ments from figure 18(b) corrected for aspect ratio and 
taper ratio by the factors given in 18(c). 



APPENDIX B 



TABLE OP AIRFOIL CHARACTERISTICS (TABLE I) 

A form of presentation of the airfoil characteristics 
has been adopted that permits all the characteristics 
necessary for the solution of problems such as those 
considered in this paper to be compactly presented. 
All such characteristics for a given airfoil section are 
presented by entries across a single line of a table. 
Characteristics are given in this form for certain well- 
known and commonly used airfoil sections in table I. 
The information presented for each section is discussed 
in the following paragraphs under subheadings corres- 
ponding to the column headings in the table. 

Airfoil: The first column of the table gives the com- 
monly used designations of the airfoil sections. 

Reference: The second column gives the reference 
to an N. A. C. A. report or technical note (R or N), in 
which additional data for the section, including the 
official table of ordinates, may be found. 

CLASSIFICATION 

Chord: The letters in this column classify the air- 
foil sections with respect to the type of their chord. 
The letter A designates a chord joining the extremities 
of the mean line, the N. A. C. A. 2412, for example; 
B designates the chord as being tangent to the lower 
surface ; and C designates an arbitrary chord from which 
the section ordinates are specified. 



PD: The letters and numbers in column PD classify 
the airfoil section with respect to the character of the 
pressure distribution about the section. The letter 
refers to the character of the additional and the accom- 
panying numbers to the character of the basic pressure 
distribution. The section of the present paper that 
discusses load distribution over an airfoil section, 
figures 6 and 7 that give the various distributions for 
the airfoil classes indicated, and the sample calcula- 
tion of the pressure distribution about the U. S. A. 
35-A section should be referred to for further details. 
The typical pressure distributions employed are based 
on Theodorsen's method (reference 7) modified to 
improve the agreement with experiment. The modi- 
fied method and some experimental results may be 
found in reference 12. No data are available for class 
A airfoils. 

SE : The character of the scale effect as affecting the 
maximum lift coefficient is indicated by the classifica- 
tion in column SE. The numbers and letters correspond 
to the designations of the typical scale-effect curves 
presented in figure 20 except that no data are available 
for class A airfoils. This information is necessary for 
determinations of stalling speeds. The Reynolds Num- 
ber corresponding approximately to the stalling speed 
is first determined. Then from the curve of figure 20 
corresponding to the designation in the SE column, the 
increment AC Lmax corresponding to this Reynolds Num- 




1,000,000 



5 6 10,000,000 6 /,0OO, 000 
fteyno/ds Number- 



10,000,000 



To obtain the section maximum lift coefficient at the desired Reynolds Number, apply to the standard-test value the increment indicated by the curve that corre- 
sponds to the scale-effect designation of the airfoil. For type A, Acj =0. 



Figure 20.— Scale-effect corrections for ci ma 
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ber is obtained. This increment is added to the stand- 
ard test value of the maximum lift coefficient given in 
table I to obtain the maximum lift coefficient to be ex- 
pected for the particular airfoil section in flight at the 
Reynolds Number corresponding to the stalling speed. 
Application of the section data to the prediction of the 
0 Lmax of tapered wings may be found in reference 2. 

This method for the prediction of maximum lift coef- 
ficients in flight is based on scale-effect tests of a number 
of related airfoils. The experimental data and a more 
complete discussion of the subject may be found in 
reference 13. 

G Lmax : Under the heading G hmax the airfoil sections 
are classified according to the character of the lift- 
curve peak. The airfoils are classified A, B, C, and D 
in accordance with behavior in the neighborhood of 
maximum lift. 

In type A the lift is more or less steady until it breaks 
suddenly to a lower value without an appreciable 
change of angle of attack. 

In type B the lift becomes so unsteady and erratic 
as to preclude the taking of measurements for a range 
of angles of attack beyond an angle referred to as that 
of maximum lift. 

In type C, before reaching the lift referred to as the 
"maximum/' the lift breaks intermittently from a rather 
definite value to another rather definite but lower 
value, and then returns to the higher value. As the 
angle of attack is increased, the breaks become more 
frequent and of longer duration. The maximum lift 
is taken as the higher value occurring at an angle of 
attack at which it is a maximum or beyond which the 
higher value can no longer be determined with con- 
fidence. 

In type D the lift is reasonably steady in the neigh- 
borhood of the maximum or any breaks occurring are 
small so that average values of the lift are measured 
throughout the range and the lift coefficients are repre- 
sented by a continuous curve in the neighborhood of 
the maximum. 

FUNDAMENTAL SECTION CHARACTERISTICS 

Effective Reynolds Number: The values in this 
column represent the values of the Reynolds Number 
at which the section characteristics should be con- 
sidered as applying to flight. The effective Reynolds 
Number is obtained from the actual test Reynolds 
Number by the application of a factor to allow for the 
effects of turbulence present in the tunnel. The tur- 
bulence factor 2.64 has been used for the variable- 
density tunnel. Comparative tests (reference 14) indi- 
cate that, at the effective Reynolds Number, maximum 
lift results from the tunnel tend to agree with those in 
flight. 

Ci : This column gives the maximum lift coefficients 

'"max a 

corrected to represent values for the airfoil sections. 

ai 0 : In this column are tabulated the angles of zero 
lift in degrees. 



a 0 : This column gives the slope (per degree) of the 
curve of lift coefficient against section angle of attack, 
that is, the lift-curve slope for a section of a wing of 
infinite span. The corresponding slopes for wings of 
finite span are found from the a 0 values by the method 
indicated in figure 12. 

Ci opt : The optimum lift coefficient, that is, the lift 
coefficient corresponding to the minimum profile-drag 
coefficient for the section, appears in this column. 
The profile-drag coefficient for the section at any lift 
coefficient may be inferred approximately from c lgpV 
c dn , and Ci by the method indicated in figure 2. 

u min max 

c dn : The values in this column give the minimum 

"min 

profile-drag coefficients. The values given, however, 
are not the ones read from the usual plots of profile- 
drag coefficient made directly from the test data. They 
are corrected for the different skin-friction coefficients 
to be expected at the effective rather than at the test 
Reynolds Number (see footnote on p. 21 of reference 
13) and for support interference. The support- 
interference correction, which gives an important 
reduction of drag for the thicker airfoils, was evaluated 
only recently and results published heretofore do not 
include the correction. Furthermore, another small 
correction is applied to these data to allow for a tip 
effect present in the tests of rectangular-tip airfoils. A 
corresponding correction has been applied to certain 
other characteristics including a 0 and the maximum lift 
coefficient; other characteristics are indirectly affected. 
A discussion of this subject may be found in reference 
13. 

c ma e : The values in this column give the pitching- 
moment coefficients referred to the aerodynamic center 
of the section rather than to the usual quarter-chord 
point. The aerodynamic center, by definition, is the 
point about which the pitching-moment coefficient is 
constant. Experimental results indicate that, by the 
use of an empirically derived aerodynamic-center posi- 
tion as suggested by Diehl, a constant pitching-moment 
coefficient c ma c may be specified for each section 
that does not depart from the measured pitching- 
moment coefficients by more than the experimental 
error, over the range of lift coefficients between zero 
lift and slightly below maximum lift. 

a. c: In these two columns the coordinates of the 
aerodynamic center ahead of and above the quarter- 
chord point are given in percentage of the chord. 

DERIVED AND ADDITIONAL CHARACTERISTICS 

Cj lc dl , : The values of this ratio are given because 

max u m { n 

the ratio has been employed as a speed-range index. 
Strictly speaking, for this purpose, values of Ci max and 
c dn should not be taken at the same value of the 

Reynolds Number; but the method has the advantage 
of simplicity and is of some value in comparing airfoil 
sections. 
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c. p. at Ci max : Values are given in this column repre- 
senting the center-of-pressure position in percentage of 
the chord behind the leading edge, or the forward end 
of the chord. The values are the measured values. 

Wing characteristics .A =6: Wing characteristics are 
given for a wing of aspect ratio 6 having the given airfoil 
section and for a modified rectangular plan form with 
rounded tips. (Tip length approximately one chord.) 
The values of m 6 represent the slope of the curve of lift 
coefficient against angle of attack expressed as changes 



in lift coefficient per radian. The values of G Dmin 
represent the minimum drag coefficients for the wings. 

Thickness. — Data are given in three columns that 
refer to the airfoil section thickness at the indicated 
representative stations. The thicknesses are measured 
along perpendiculars to the chord and are expressed in 
percentages of the chord. 

Camber. — The camber expressed in percentage of the 
chord is represented by giving the maximum displace- 
ment of the mean line from the straight line joining its 
extremities. 
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SYMBOLS 
BASIC CONSIDERATIONS 

S, wing area. 
L, lift. 

g_, dynamic pressure (l/2pF 2 ). 
c, chord. 
b, span. 

^ L== qh' wm & ^ coefficient. 
dL 

Ci= — j~> section lift coefficient. 
gcdy 

c l() , section lift coefficient acting perpendicular to 
local relative wind. 

Subscripts: 

al, refers to additional part of load distribution 
for C L =1. 

a, refers to additional part of load distribution 

for any C L . 

b, refers to basic part of load distribution for 

C L =0. 

y, distance along lateral airplane axis. 
general procedure 

Monoplanes 

x, distance along longitudinal airplane axis. 
z, distance along normal airplane axis. 
Xa.c, 2a. c, x and z coordinates of wing aerodynamic 
center. 

L a , additional load parameter, c lal c^- 

J.. j (3 j J.- and z coordinates with respect to a system 
of axes originating at the wing aerodynamic 
center. 

L b) basic load parameter, c ib — ^- 

X, Z, components of air force in the x and z direc- 
tions. 

c x , c z , section force coefficients. 

M, wing pitching moment. 
Ma.c.i wing pitching moment about wing aero- 
dynamic center. 
M s , part of wing moment due to section pitching 
moments about their aerodynamic centers. 
Mx, M z , parts of wing moment due to X force and Z 
force. 

M T , wing pitching moment about torsional axis. 
Ztj Xj', distances of the torsional axis below the chord 
plane through the section aerodynamic 

22 



center and behind the beam plane through 

the section aerodynamic center. 
c ma c , section pitching-moment coefficient about 

section aerodynamic center. 
c dQ , section profile-drag coefficient acting parallel 

to local relative wind. 
a 0 , section lift-curve slope (per degree). 
ai 0 , section angle of attack for zero lift. 
i, incidence of section chord with respect to x 

axis. 

c„ tip chord (for rounded tips, c t is the fictitious 
chord obtained by extending leading and 
trailing edges to the extreme tip). 

c s , chord at center of wing or plane of symmetry. 

A, aspect ratio, b 2 (S. 

e, aerodynamic twist, assumed linear, and meas- 
ured as the angle between the zero lift 
directions of the center and tip sections, 
positive for washin. 
C L ' ', wing lift coefficient for steady-flight condition 
preceding accelerated-flight condition, 
w, slope of wing lift curve (per radian). 
m 6 , slope of lift curve for nontapered wing with 
rounded tips and aspect ratio 6 (per radian). 
c<7„ , minimum section profile-drag coefficient. 

c , , optimum section lift coefficient, lift coefficient 

corresponding to c rfo ^. 

V, flight velocity or air speed. 
U, velocity of gust. 

a, slope of wing lift curve (per degree). 
/, plan-form factor. 

B, beam component of force. 

C, chord component of force. 

c b , c c , section coefficients of beam and chord com- 
ponents. 

Uq 

(f>=i b —i c 

i b , incidence of section chord with respect to the 

perpendicular to the beam direction. 

i c , incidence of section chord with respect to the 

chord-truss direction. 
C N , wing normal-force coefficient. 
c n , section normal-force coefficient. 
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Biplanes 

x, z, distances defining the locus of the aerodynamic 
centers of the biplane cellule. 
K h K 2 , etc., Diehl's biplane lift functions (references 5 
and 6). 
G, gap. 

t, thickness of wing. 

s, stagger, distance between aerodynamic cen- 
ters of upper and lower wings measured 
parallel to x axis. 
s 0 , stagger, distance between M-chord points of 
upper and lower wings measured parallel to 
the zero-lift direction. 
M Y y, net biplane pitching moment about arbitrary 
y axis. 

^ Cm °(g)» increment in section moment coefficient due 

to biplane parameter tjG. 
Su', portion of upper wing area to which the tjG 

moment correction applies. 
S L ', portion of lower wing area to which the tjG 

moment correction applies. 
c v ' , average chord of the portion of the upper wing 

corresponding to Su'. 
c r J , average chord of the portion of the lower wing 

corresponding to St! . 

Subscripts: 

U, refers to upper wing. 
L, refers to lower wing. 
B, refers to biplane. 

LOAD DISTRIBUTION OVER AIRFOIL SECTION 

P, normal-pressure coefficient, p/g. 

f, the pressure difference across the wing section 

at any station along the chord. 
P , value of P for the pressure distribution at 
zero lift. 

P a , value of P for the additional part of the pres- 
sure distribution when the additional sec- 
tion lift coefficient is 1. 

x 

PaV - ? -~P a .c, components of P a . (See fig. 6 and 
c 

equation (24).) 

-^^j distance in terms of chord, of section aero- 
c 

dynamic center forward of quarter-chord 
point. (See table I.) 
P b , value of P for the basic part of the pressure 
distribution. 

—c Ma c P bm , ~ P bc , components of P b . (See fig. 7.) 
c 

— i camber ratio; distance, in terms of chord, of 
the minimum height of the section mean line 
above the straight line joining its extremi- 
ties. (See table I.) 

c nb , section normal-force coefficient corresponding 
to basic pressure distribution. 

~c ma c _(c n ) bm , ~(c n )bc, components of c„ b . (See fig. 7.) 



SAMPLE CALCULATION 

W, weight of airplane. 
F t , force on horizontal tail surfaces. 
n u applied wing load factor, basic design Condi- 
tion I (a), reference 1. 
s, effective wing loading, reference 1. 
(Jl, dynamic pressure corresponding to design 
high speed. 

C Nv wing normal-force coefficient corresponding 
to n v 

a s , wing angle of attack based on chord of central 
section. 

(aio) s > angle of zero lift of the central section. 

J, parameter for determining angle of zero lift of 
twisted wing. (See fig. 11.) 
Cd 0 ' t section profile-drag coefficient for steady- 
flight condition. 
q', dynamic pressure for steady-flight condition. 

Biplane 

I, distance from y axis to c.p. of horizontal 
surfaces. 

F x , F 2 , factors for reducing wing lift coefficient to 
allow for tip increment. (See fig. 17.) 
k, Munk's span factor. 

REFERENCES 

1. Bureau of Air Commerce, Dept. Commerce: Airworthiness 

Requirements for Aircraft. Aero. Bull. No. 7-A, U. S. 
Dept. Commerce, 1934. 

2. Anderson, Raymond F.: Determination of the Character- 

istics of Tapered Wings. T. R. No. 572, N. A. C. A., 1936. 

3. Pearson, H. A.: Span Load Distribution for Tapered Wings 

with Partial-Span Flaps. T. R. No. 585, N. A. C. A., 
1937. 

4. Pearson, H. A.: Empirical Corrections to the Span Load 

Distribution at the Tip. T. N. No. 606, N. A. C. A., 1937. 

5. Diehl, Walter S.: Relative Loading on Biplane Wings. 

T. R. No. 458, N. A. C. A., 1933. 

6. Diehl, Walter S.: Relative Loading on Biplane Wings of 

Unequal Chords. T. R. No. 501, N. A. C. A., 1934. 

7. Theodorsen, Theodore: Theory of Wing Sections of Arbi- 

trary Shape. T. R. No. 411, N. A. C. A., 1931. 

8. Garrick, I. E.: Determination of the Theoretical Pressure 

Distribution for Twenty Airfoils. T. R. No. 465, N. A. 
C. A., 1933. 

9. Jacobs, Eastman N., Ward, Kenneth E., and Pinkerton, 

Robert M.: The Characteristics of 78 Related Airfoil 
Sections from Tests in the Variable-Density Wind Tunnel. 
T. R. No. 460, N. A. C. A., 1933. 

10. Bxireau of Air Commerce, Dept. Commerce: Design Infor- 

mation for Aircraft. Aero. Bull. No. 26, sec. 7 (D), 
U. S. Dept. Commerce, 1934. 

11. Diehl, Walter Stuart: Engineering Aerodynamics. The 

Ronald Press Co., 1936, p. 34. 

12. Pinkerton, Robert M.: Calculated and Measured Pressure 

Distributions over the Midspan Section of the N. A. C. A. 
4412 Airfoil. T. R. No. 563, N. A. C. A., 1936. 

13. Jacobs, Eastman N., and Sherman, Albert: Airfoil Section 

Characteristics as Affected by Variations of the Reynolds 
Number. T. R. No. 586, N. A. C. A., 1937. 

14. Jacobs, Eastman N., and Clay, William C: Characteristics 

of the N. A. C. A. 23012 Airfoil from Tests in the Full- 
Scale and Variable-Density Tunnels. T. R. No. 530, 
N. A. C. A., 1935. 



24 



REPORT NO. 631 — NATIONAL ADVISORY 



COMMITTEE FOR AERONAUTICS 



Pill 



II 



J\ i 



gggS2 8SfegsaS2£SSSSSSSS8SSS88 



88S8S&8g8Se8S3fiiS8RS889S8!8S:( 



lIlIIHIllllllillllllllllll 



T|i -* ^ ^ - 



Si 



.'I 



l§S§ig8811s21 IIsl!lSil§li 
f r f r i r i' i' ' i' r f r°° r r r i 1 1 i 1 1 i i 



illlllllililllllllllllllll 



I|gl8l§lllglll§l§llll§l§l8l 



fffffffT ''tttf 



I I I I I I I I I 



aSododt-ajodadenododododododododwadw^ 



222222222 2222222 22222222£!SS; 
OQWOfiPOQO<1DHO^OOMOQWO«!(riOQWfa 



i 



llilliilllllii 



pq : i : : i T ; i : i~7 i :JJJ JjJ 
o2 i i j Jlllllllllllll 

ddf*-«.«j«jdddddddddddddd 



THE PREDICTION OF AIR FORCES AND THEIR DISTRIBUTION ON WINGS 



25 



Si£3ssg!g§38§§ 

ggggggggggggg 
S53S18§S§SSSS 
J5S8SgS332SS3§ 

i-ir-ir-ir-ir-ii— ii-ir-ii— ii— ir-4f-ir-< 

^^^^^^^^^^^^^ 

mmmmwm 

93§38SgSSggi*g 



A 



sSiiiiigisgiii 

d 



ig§iigi£ggg§§ 



§§ii§iglggg£g 



d 



s§S8§S§8§1SS1 



d 



d 



d 



d 



Ilslllllll™ 



d 



1 



gg§3g§3i533g§g 

J^J^Wr-i^r-iwr-ir-ir-ir-ir-; 



rn ^ r4 >4 r-i .-i -i r-i ~i .-i r-i r-i r-i 



£SgggS§gSS2§S 



r-ir-ir-ir-ir-ir-ir-ir-i-i-iwr-ir-i 



HHi-ii-iHi-irtrt'rtHrHWi-l 



_i W r-i W r-i rJ r-< r-i r-i ~ -i r-i 



r^r-ir-ir-ir-ir-i^r-ir-ir-ir-ir-ir-i 



-ir-ir-ir-ir-ir-ir-ir^r-i-ir-ir-; 



H H t-l i-J H f-j r-i H H H r-I i-H 



jjjjjj 



S!§3B§S2§||iiS 



r-i ^ ^ ^ ^ ,_; ^ ^ ^ ^ ^ 



SSSSSSSSfllli 

rJ-ir-i^r-i-irJr-ir-ir-ir-ir-;^ 



_i _j „ ^ ^ „• _j r-i rJ r-i r-i 



^ir-i^.-i^^^wr-ir-ir-ir-;-; 



5S5S3:!g3S53aS 

„' r-i -i -1 -i rJ r-i r-i « r-i r-i r-i r-i 



SSSSS§222222§ 

H r-i rl r-i r-i r-i r-i -i r-i r-i 



gggSSSSggSggg 

r-i r-i r-i _i r-l r-i _i r-i rJ r-i 



r-i r-i r4 -i r-i r-i r-i r-i r-i ri ri r-i r-i 



i i i ; i i 

i i i i : ; 
:::::: 



i 



gg2§g§g§3£§§§ 



SgSSSSSSSSSSK 
d 



BSIIlggllllll 



t 

! 



11II11IISSIS1 

-iw^r-irrr-ir-ir^r-Jr-ir-i-ir-i 



sllll§llllll§ 



slllglllllsss 

HriririririrlHririHriri 



Issslsssllill 

Hriririri-i-riririri ' ' 



isiiiiiiiiiii 

^ ^; ^ ,h „; 



lilillilgllll 



llllliSllllll 



§igg§§§££§3§§ 



immmmmm 
© 



iiiiiiii§iiii 



: : 



22 = 



§£§S™£§111§ 



ggggggESSgggg 



ggggSgs£sl88S 



ggggg»g£s§ii§ 



5gg£sgsss§§g§ 



gg£g£§558§ll§ 



Sgg§£g£gg§§l§ 



ggggggsssggll 



§h 



rrr 



1 1 

§lll§gg§§l°§l 

o I I 1 I I o 

§§!§§§ §§!§gg 

o I I I I I I I I d 

SSS 888233222 

a ^ to MtooooooOH'Mw 

I I I I I I I I I II I o 

§§§§§§§§§§§§§ 

I I I I I I I I I I I I o 
I I I I I I I I I I I I o 

S8S88S2222222 

f I I I I I I I I J I I o 

§g§§l§g§§§§§§ 

I [ I I [ I I I I I I I o 
I I I I I I I I I I I I o 

I I I I I I I I I I I 1 1 

I I I I I I I I I I I I o 

§§!§§§§§§§§!§ 
g§§§§§§§g§§|j 

Sg'gg'g§§§883il 

'5g'g'ggg§g§§l§l 

ggggSg'ggg'ggll 
gggggSgggglil 
gg'ggggg'gggggl 

sgggggggggggi 
ggggg§gggg§§l 

Sgggggggggggl 



p 



I I I I I I I I I I I I o 

gggiggggggjl 

I I I I I I I I I I I I o 

I I I I I I J I I I I J <L 

I I I I I I I I I I I [ o 

Sgggggjggs'gg'sg 

I I I I I I I I I I I I o 

I I I I I I I I I I I I o 

SgggggiSSgggg 

I I I I I I I I I I 1 1 o 

g'gggggggS&sgg 

I ! I 1 I I I I I I I I o 

ggggggggsggsg 

I I I I I I I I I I 1 I o 
I I I I I I f I I I I I o 

§gjgS§g§ §agggS 

I I II I I I I I I I I o 

§8§§g §§§§gl§l 

Sggg'5'ggg'5§§§l 
p 

"ggggggggggggl 
gg ggggggggg ll 

ggsggggg'sk'skl 
gggg ggg'g§g g§l 
gggg'ggggg'g'ggl 
ggggggggggggl 

ggggggg'aggggl 

ggggggggggggi 



§1* 

ii 



I I I I I I I I I I I I o 

gg§§§Sggggggs 

I I I I I I I I f I I I o 
I I I I I I I I I I I I o 

5§ggggggggg§B 

I I I I I I I I I I I I o 

Si§glgggggg§gg 

I I I I I I I I I I I I o 

gggggg'ggggggg 

I I M I I I I I I I I o 

igsgggggggggg 

I I I I I I I I I II I o 

gggggggggggg's 

I I ! I I I I I I I I I © 

ggggssggggggg 

I I I I I I I I I I I I o 

I I I I [I I I I I I I o 

gggggggggg ggg 

I I I I I I I I I I I I o 

ggggggggggggg 
gggggggggg ggj 

gggg'a'g'g'g'ggggl 

g§§gg§gg'ggg§i 

gggggg's'g'ggggl 
ggggggggggggl 
gggggg'sSggggl 

ggggggggggggl 
ggggggggggggl 
ggggggg'sggggl 

ggggggggggggi 



V 



f 



W 

CO 

l-H 

O 

CO 

3 



o a 

S H 

u o 

tri M 

U h 

„ a 



g§g§s§§§§§§§§ 



g8g§S8§§ll§i§ 



gggggs§gg§S§§ 



gggggg§gg§ii§ 



ggggggsgll§g§ 



ggs5gg£s§g§§§ 



■Sh 
ii 

p 



Sh 
p 
s 



Sh 
p 



Sh 
p 



o 
> 



soiiLnvNonav ho..* aaxxiKwoo ahosiaciv tvnoixvn— is9 on x'aoaara 



9S 



CO 

o 

O 

o 

t— 1 
Eh 
P 

m 

H 

£h 

GO 
I— I 

ft 
P3 

p 

H 

n 

02 

w 

O 
P4 
O 
P 

< 

o 
o 

t— I 

H 

O 
i — i 

ft 

H 

Ph 

H 
K 
H 



O 
H 

O 

ft 
<J 

P 

PK 
^O 

Oq 

ft<^ 

m 9 

HO 
D# 

o< 

£P 
ft^ 
ftp 

<io 

P£ 
wo 

p<: 

Op; 

^P 
pft 

3* 

!> 
i— i 

H 

<J 
H 



2 


as 
1 


-0. 0366 
-. 0175 
-. 0035 
.0070 
— . 0035 
-. 0052 
0 

.0035 






T - 


MM 




1 


0. 999 
1.000 
1.000 

1. 000 
1. 000 
1. 000 
L 000 
1.000 




S3 


u ST- 


J 




n 




3 




, 1 1 II 






n 


0.416 
.248 
. 128 
.043 
.077 
.063 
.033 
.017 


OS 


jo" 


0. 0157 
.0159 
.0157 
.0147 
.0124 
. 0117 
.0115 
.0115 




s 




l l 1 1 


X 


T | 


TTTTTTTT 




2 






I s - 


•1 


111111111 


2 




0. 942 

1. 000 
1.035 
1. 03fi 

.945 
.911 
.924 
.942 


S 


•3 


iiiliiii 






i 




1 


-0. 0006 
-.0003 
-. 0001 
.0001 

0 

—.0001 
0 
0 


« 


J 

s 


IIIIIIIII 
© 






0. 944 
1.000 
1.035 
1.030 
.945 
.911 
.924 
.942 






1 
« 


0. 943 
1.000 
1.035 

1. 036 
.945 
.911 
.924 
.942 


CO 








■s 

J? 

T 


linn i 

° i © i 






1. 301 
1.240 
1. 142 
1.000 
.775 
.583 
.436 
.320 

0 


2 


1 

■? 


0.0157 
.0159 
.0157 
. 0147 
.0124 
.0117 
.01 15 
. 0115 




Station from 
center line 
y 

b/2 






Station 
from cen- 
ter lino 

b/2 


o 



P 

O 

2p 

Eh 

33 

P H 

O 

< 



<1 
> 



3 



;So32o 



Iiiliiii 



d 



J 




J 


d 




Illllilll 




Illllilll 




o 





fl 


IHsslls 
° 1 1 11 1 






Tftt T7 1 


00 


n 


Iliiilii 

d 




S 


-0. Ill 
-. 104 
-. 097 
-.090 

— . 083 

-.078 
-. 077 




1 


2011 
1925 
1891 
1822 
1771 
1754 
1719 
1719 


Si 




^ io o x « m ^ « 


<= 1 II 1 1 1 1 


o 




if 


!§§!!!!!! 

o 


8 




-0. 08 
-.02 
.02 
.02 
. 01 

!oo 

.00 
. 00 


2 


4 


ri32222~f 

i i i i i i 


S3 


g 


1 1 






liiiiiii 


CN 




,H |' |" ' ~J r-«' IN 


O0 




TTTTTTTT 






° Mill 


r~ 






11111111 


2 




d i 


o 




° 1 1 1 1 1 






liiiiiii 




fl 


liiiiiii 


° Mill 


d 






iliiilii 
= r ■ 






0. 0400 
,0312 
. 0040 
-.0231 
-.0427 
-. 0530 
-. 0582 
-. 0540 






-0.0027 
-. 0024 
-. 0023 
-. 0021 
-. 0019 
— . 0018 
-. 0017 
-. 0017 


CO 


1 


00 as o M N ?i 






T 1 1 


3 


x ^ 

-1 


Iliiilii 

° ,11m 




Station 
from 

center 
line 

b/2 


o 




Station 
from 

center 
line 
V 

b/2 


!M « 00 CS S £ 

O 



O 

Oh 

C 

o 

< 
pq 
Q 

G 
Q 

£ 

H 
EE 

o 
o 

I — I 

H 



> 



00 


sin 0 /> 


mmmm 


r- 


1 


0. 984 
. 979 
.977 
. 970 
! 981 
.983 
.982 
.981 


s 


= 


Ottilia 




5 


d 




li 

3 


o «s 00 cd i-i o 
co co ci cA cm' m <Ni 




3 


■5T >r.' IO CO C*i CC CO 




-a bo 


o -.' n c4 d d -! 




■fi 

2. 


CO CO c<i c4 c4 ri 


2 


f| 


oooooooo 


cs 


r | 


od t-: d d >rf o ^ ■« 
MMMM 


00 




Wr-CCOCiONO 

°S222S5S3 


t_ 




0. 0120 
.0110 
.0110 
.0104 
. 0093 
. 0090 
. 0087 
. 0085 




2° 


liiiiiii 

d 


m 




-0. Ill 
-. 104 
-.097 
-. 090 
-. 083 
-. 079 
-. 078 
-.077 


■* 


<? 




CO 


i 


1.702 
1.803 
1.865 
1. 860 
1. 703 
1. 642 
1. 606 
1. 698 


CS 

I 




ummm 
° M 1 1 1 


! 


Station 
from 

center 
line 

b/2 


c^oooogl 



05 



O 

2; 



2 
o 

o 

I— ! 

E- 

i — ■ 

i — i 
P 

K 
K 
E- 

P 



C 



ft 

o 
c 

H 

o 

I— i 

p 
a 



Eh 

H 

o 



> 

g 

CH 



ft," 




sf.! 

c J 

i „ 15 
ft, tt, 


° 11111° 


ft," 


o to «i to ui ui «i « -! ■ p4 <M' <o »o o 
Mill 


■3 
ft< 


S5£S53$|g3pS8g3SSi = S 




o^ci^oooooooooggg 



o 



a 

o 
o 

1—1 

H 
<1 
(-1 

o 

<; 
o 



2 


-sSsgsstfasSSp^- i 


ft, 




ft," 
s 






SSSaSsSpiSSgSgSSSpp 
MM 


t- 

T 


S§£g£3SS3gSc$pt§3 

°7 7TTTT7 1 ' '° 


aT 


O pH i-H O 


ft. 




s 
ft. 


ooooooooooooooooo 


i 

"j 
ft, 
1 


O r-<"rH o 


4 


S38gS8ggSgiS8S§i?22 


Station 

100-? 
c 


So a 

o^««r»:o«oooogoogog . 

1 



Rib location in percent semispan 


Si 


Aspect ratio 4 




Aspect ratio 0 




Aspect ratio 8 


Sgt2ggg«Sp?^^ 







pM 


i§ 






3§§SES383i§ 

PHM 


g 








8 


iisiiiisiig 




IsSlllllIIs 




pH p-i -M ph" ph - 


piripipip'pi 


g 


plSllgglllll 

^ w P-i r-l „, w ph' _i 


sSillllllll 

_• pH ph pH rH - Hri pH r-Jr-i 


SS§S§§1ISS8 

-irtrtW PH* P- PP PH P-i pI 


S 


lllllllls^ 

_• „• PH ^ PH- p_j _i ^ PH 


lllllllllss 

_! ^ h H ph" P-i pH pH 


iliiiiislis 

j ph ph - p-; p,' ph- p,- 




Ss3S§gS83SS 

O • ' ■ p-' pH pH r-i pH pH pH 


isigiiisisfp 


llllgglg§§s 

o ' ' ■-•pi-ipiHH'p; 


a 












O pipipipi 


Taper 
ratio 








o ph 


o p.. - 


o i — 1 



o 

I— I 

o 

Ph" 

< 

M 
O 

Ph 



o 
o 

P3 

o 

CO 

i— i 
> 

n 
< 

i— ! 
< 

O 

HH 

<5 



O 

o 

Ph 

H 



c 
o 

hH 

Eh 
P 
pq 



ft !i 

O ^ 

H a 

H 2 

1 

O "3 

a s 

Eh 





U3 

s 


rfrt^H^H^-. I <N (N CO OC 

f 1 1 1 1 1 1 ! 1 1 1 












^ i i i i i i i i i i 




§ 






=j> 1 1 1 1 M M II 






f 1 1 1 1 1 1 1 1 1 1 


i 




°i 1 1 u 1 1 1 1 1 


location in p 


§ 


!§§§§8S§SSS 
°lllllillll 




s 


giiiigiiiii 
° 1 1 i 






0.043 
.040 
.037 
.034 
.030 
.026 
.023 
.019 
.010 
. 012 
.008 






0. 165 
. 158 
. 154 
. 149 
. 145 
. 140 
. 134 
. 126 
. 117 
. 103 
.083 




o 


o 


If 


c — 



Ph 

hH 

pq 

Ph 
o 

a 

Z5 



Pi 

Cm 
Ph 

B 

o 

o I 

t— I ^3 

H S 
p 1 

PQ * 

hH S 

P3 | 
§ 1 

q & 
h 

o a 

P3 -S 

o I 

H 

a 

H 

o 
o 

Eh 
< 

Ph 

C 

o 

R 

H 

pq 
< 





+ 

r s 


0. 0075 
.0075 
.0075 
.0074 
.0073 
. 0073 
. 0073 
.0072 
.0071 


El 


1 


S88SSS1SI 

1 Ml 




i 


iillllll ! 




1 


3ggg333Sa 




° i i i i i i u 


o 


1 




§ 


< 


o o o 

OOH 

C O O O O O |" |' [' 


C5 


i 


0. 2111 
.2097 
. 2045 
. 1906 
. 1605 
. 1395 
. 1300 
. 1340 
. 1449 


2 


f 
< 


o ' ' ■ o o o o o 




s 

o 




2 


r 


IllSIIIIs 
^Iillllll 




S 


12.2 
12. 1 
11.8 
11.0 
9.2 
8.0 
7.5 
7.7 
8.3 




s 




eo 






7777 I l l 1 l 
MM 


2 


1 


gllllllll 

O 




+ 

S 


1. 357 
1. 349 
1.319 
1.240 
1.070 
.952 
.899 
.921 
.983 


2 


OS 

■1 


seisms 

7^77 Mill 




1 


o o o o c o 


CO 


ji 
V 


§133§§£g§ 




I 
+ 

2 








1. 137 
1. 130 
1. 105 
1.039 
.897 
.798 
.735 
.638 
.480 


2 


1 


oooooooog 




Distance 
from center 
line (ft.) 






Distance 
from 
center 

line (ft.) 





o 
co 




tv ^ 

A Diameter 
P, ' Geometric pitch 
1 2>/A Pitch ratio 
V '> : ' Mow: velocity 
V», Slipstream velocity 



4. PROPELLER SYMBOLS 

Ilfl 



Power, absolute coefficient CV==— 

■ - Speed-power eoefficient-^i/^ 
^ Efficiency 



• Thrust > absolute coefficient <?r=^ ■ ■ Resolutions per second, r.p.s. . ^ 

G Torque, absolute coefficient C«~J^ ' • ffe ° tn 6 hellX a ^ Ie==t ^ 

U D _ 7fifti , ' - - 6 ' NU MERICAL RELATIONS 

1 Hp 76.04 kg-m/s«550 ft-lb./sec. 

1 metric horsepowers 10132 hp ih.=0.453G kg. 

1 m.p.h.==0.4470 m.p.s. 1 %=2.2046 lb. 

1 " 1 mi - = 1^09.35 m=5,280 ft 

1 m=3.2808 ft. 



1 m.p.s. =2.2369 m.pji. 
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